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-Th«  Any  Multiple  Systems  Evaluation  Program  (MSEP)  Is  a 
oosprehenslve  program  developing  general  analytic  techniques 
for  the  prediction  of  ‘hlgh-electromagnetlc-pulse 
vulnerability  and  hardening  technology  and  for  the 
Application  of  these  techniques  to  a list  of  critical 
systems.  Ihe  analytic  techniques  have  been  verified  for  a 
large  class  of  tactical  systems.  The  hardening  techniques 
have  been  applied  to  specific  systems  and  are  now  resulting 
in  product  improvement  programs  leading  to  hardened 
equipment  in  the  field. 
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1.  XMTIKMDDCTZOM 

Thi*  «ffort  wu  writton  undar  the  sponsorship  of  the  ttoltiple 
SystSM  Bvslustlon  Prograsi  (MSEP)  iriilch  has  as  its  main  objective  to 
harden  hzay  tactical  systsa«  to  the  exoatno^heric  electrcaiagnetic-pulse 
CMP)  threat.  Along  with  this  najor  objective,  the  MSEP  is  also  tasked 
with  the  ala  to  develop  experiaental  and  analytical  evaluation 
taohniques  that  are  applicable  to  all  systems'  problesis.  This  work 
illuadnatas  several  numerical  techniques  that  can  be  used  to  solve 
ordinary  differential  equations  that  night  arise  in  making  EMP 
vulnerability  assassnenta  for  tactical  Amy  systems.  Therefore,  the 
objective  of  providing  analytic  methods  is  satisfied. 

A study  was  undertaken  to  expand  «qpon  the  results  obtained  in  a 
paper  by  Myatt^  and  to  investigate  and  compare  two  numerical  techniques 
in  solving  ordinary  differential  aquations  that  occur  in  an  EMP  problem. 
Zt  was  conducted  with  the  purpose  of  refining  the  previous 
calculations.^  The  numerical  tedmlques  employed  are  fairly  standard 
methods,  but  whereas  Nyatt  was  concerned  only  with  approximate 
solutiona,  this  work  deals  mainly  with  the  qpeclfic  numerical  techniques 
uaed  In  the  solution  of  the  problem.  An  analysis  of  truncation  errors 
was  employed  to  deal  with  the  accuracy  of  the  solutions  that  better 
solidifies  the  results.  The  methods  utilised  were  the 
prsdlcbor-oorrector  technique  known  as  HsMlng's  method  and  an 
emonentlal  differencing  method  that  was  developed  by  Pope.^  This 
effort  does  not  attempt  to  esqplaln  the  physics  of  the  problem  or  the 
derivation  of  the  ordinary  differential  aquations.  This  material  is 
fully  detailed  by  Nyatt 

Zn  the  following  sections,  soam  general  comments  are  made  concerning 
Hunge-Ktttta  solutions  and  predictor-corrector  solutions  of  ordinary 
differential  equations.  The  two  techniques  employed  in  the  solution  of 
the  pcoblma  are  described  fully,  along  with  a discussion  and  ccsqparison 
of  the  results.  At  the  end,  some  conclusions  are  made  concerning  the 
results  and  the  numerical  methods  that  were  used. 


^N.  t,  Npatt,  Zhteraal  IMP  Strength  and  Tima  Dmpaodmacm  for  an 
ZdmmliMmd  Pxoblm,  Jtsport  1994,  U.S.  Armg  Mebility  Bqulpmmnt  JCsssarch 
and  bsvmlopasnt  Cmntmr,  rort  Bmlvoir,  Ml  /February  1971), 

^Dmrid  A,  Fops,  An  Ixpomentiai  Matbod  of  mmarioal  Zutagration  of 
Ordinary  Diffarantial  Equationa,  CoauNinicationa  of  the  ACM  6^,  so.  8 
/August  1963),  491-493. 
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2.  GENERAL  PHYSICAL  PROBLEM 

Th«  problaa  that  was  dealt  with  tms  the  datarBlnation  of  the  tiM 
history  of  the  electric-field  strength  between  two  parallel # infinite, 
alwinusi  plates  caused  by  Compton  electrons  generated  by  a transient, 
high-intensity,  gamma-ray  flux.  The  ordinary  differential  equations 
that  are  solved  are  a result  of  the  application  of  Poisson's  equation  to 
the  free  charge  that  is  the  spatially  distributed  current  of  Ccapton 
electrons.  Since  the  conductivity  of  the  air  is  significant,  Poisson's 
equation  is  modified  by  Ohm's  law,  tdiich  relates  the  air  conductivity, 
electric-field  strength,  and  conduction  current.  In  the  case  examined 
here,  methane  was  used  instead  of  air,  and  it  was  found  by  Hyatt^  that 
methane  reduced  the  peak  electric-field  strength.  The  resulting 
ordinary  differential  equations  to  be  solved  are 

. at » (oj.'Sd)  - , t<o  u) 

CLZ  t m t 


^^dt^  + f* **-  N^(t)E(t)  - — , t > 0 , 


where 

E(t)  is  the  electric-field  strength, 

q is  the  electronic  charge, 

u is  the  electron  nobility, 

c is  the  permittivity  of  free  space, 

N^Ct)  is  the  free  electron  density, 

r is  the  model  parameter  for  gamma-flux  rate  history, 

J is  the  Cosupton  electron  current  density, 
c 


*ir.  r.  ttyatt,  IntexnaJ  BMP  Etreagth  and  Time  Dmpendencm  fox  an 

Idomliamd  fxoblem,  ttmport  1994,  V.S.  Jkxmg  Mobility  Eguipsmnt  Mesesreh 
sad  Osvnlopnsne  Center,  fort  Belvoir,  n (febmexy  1971} . 
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3.  RDNGB-KDTTA  AMD  PREDICTOR-CORRECTOR  METHODS 


Mhen  attMptlng  to  solve  ordinary  differential  equations 
maerically,  two  well-known  methods » the  Runge-Kutta  technique  zuid  the 
predictor-corrector  technique^  are  often  en^loyed  to  find  the  solution. 
Depending  on  various  constraints  such  as  time#  computing  funds 
available,  accuracy,  and  steUaillty,  there  are  accompwying  advantages 
and  disadvantages  to  both  methods.  Thus,  before  examining  HaimBing's 
method  in  the  solution  of  equations  (1)  and  (2) , a few  general  remarks 
concerning  fourth-order  Runge-Kutta  solutions  and  fourth-order 
predictor-corrector  methods  are  made  to  give  some  backgrotmd  and 
rationale  for  the  technique  utilized. 

It  is  generally  recognized  that  predictor-corrector  methods  are  more 
difficult  to  code  than  Runge-Kutta  methods.  However,  although 
Runge-Kutta  techniques  are  more  straightforward,  predictor-corrector 
methods  provide  a much  easier  analysis  and  examination  of  errors  and  are 
generally  much  faster.  For  exanqple,  the  evaluation  of  f(x,y)  (l.e., 
dy/dx  " f(x,y))  is  usually  the  most  time-consuming  part  of  solving 
differential  aquations,  and  fourth-order  Runge-Kutta  methods  require 
four  evaluations  of  f(x,y)  per  step,  while  fourth-order 
predictor-corrector  SMthods  require  only  two  evaluations  of  f (x,y)  per 
step.  This  means  that  the  fourth-order  predictor-corrector  methods  are 
generally  nearly  twice  as  fast  as  fourth-order  Runge-Kutta  techxd.ques. 
Thus,  the  evaluation  of  errors  and  the  speed  of  cooptation  axe  two 
ooa^elling  reasons  to  use  predictor-corrector  methods  instead  of 
Runge-Kutta  methods.  However,  since  predictor-corrector  methods  axe  not 
self-starting  and  Runge-Kutta  techniques  do  have  the  self-starting 
capability,  Runge-Kutta  methods  are  quite  useful  in  generating  starting 
values  for  the  solution  and  may  be  used  to  change  the  interval  between 
steps  when  desired.  This  usefulness  makes  Runge-Kutta  omthods  an 
indispensable  tool  in  using  predictor-corrector  techniques.  Therefore, 
a cosd>inatlon  of  these  two  swthods,  (1)  the  Runge-Kutta  to  determine 
starting  values  and  change  the  per-step  interval  and  (2)  the 
predictor-corrector  to  actually  solve  the  differential  equation  and 
analyse  the  errors  involved,  gives  a technique  that  is  highly  desirable 
in  coopting  the  solution  of  ordinary  differential  equations. 


4.  RAMIZIIG'8  NETHCH) 

Since  HasBlng's  method  is  not  self-starting,  the  Runge-Kutta 
technique  was  used  to  start  the  solution.  This  involves  mitlng  the 
differential  equation  as 

g - t(x.r) 
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and  than  finding 


k2  - + ytl)h  , 

k3  « f^Xj^  + + %k2)h  , 

k4  - f^Xj^  + h,  y^  + ka)h  , 


where  h is  the  desired  increment  step  size.  Once  these  vzdues  are 
computed  for  a particular  increment#  ^y  is  calculated  by 

Ay  - i(ki  + 2k2  + 2k3  + k^)  , 

and  then 

- ^i  + Ay  • 


This  process  is  continued  until  the  desired  number  of  starting  values  is 
obtained.  In  Haiming's  method,  three  values  from  the  Runge-Kutta 
technique  are  used  to  start  the  solution  process. 


Once  the  starting  values  have  bean  detendned  from  the  Runge-Kutta 
technique,  the  general  procedure  for  Hamming's  method  is  used,  as  is 
dona  by  most  predictor-corrector  techniques.  This  procedure  is 
illustrated  by  the  flow  chart  in  figure  1.  The  q^ific  formulas  used 
in  Hamming's  method  and  shown  in  figure  1 are  as  follows: 

y“;  - - yJ.j  . 2yj^)  , 


Oorzeetort 


^1+1 


Truncation  error:  T^  a, 

where  h is  the  interval  step  else. 


ly  . ytO) 

I'l+l  ^i+1 


]' 
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In  applying  Haoalng'a  nethod  or  any  pradlctor-corractor  method,  one 
must  be  cautious  about  the  number  of  iterations  of  the  corrector.  To 
let  the  iteration  of  the  corrector  continue  for  any  length  of  time  would 
defeat  one  of  the  main  reasons  to  use  predictor-corrector  methods,  qpeed 
of  computation.  It  is  an  advantage  to  know  idien  to  stop  iterating  the 
corrector.  Therefore,  it  is  reasonable  to  require  that  any  error 
by  not  iterating  the  corrector  to  convergence  be  souill  when  coo^Mured 
with  truncation  and  roundoff  errors.  Generally,  truncation  errors 
affect  the  solution  of  differential  equations  and  cause  instabilities 
more  than  roundoff  errors;  hence,  in  predictor-corrector  techniques, 
only  tmncation  error  is  considered,  with  the  above  point  in  mind,  the 
following  procedure  was  used  to  estimate  the  error  incurred  by  not 
iterating  the  corrector  to  convergence.  After  each  corrector  iteration. 


.tl) 


- y 


Ci"l)  • 


was  computed  and  coapared  with  a convergence  factor,  c (e  in  fig.  1) . 
Bie  value  e was  ctipsen  so  that  if  < c,  then  terminating  the 
iteration  with  ' and  y(i+U  would  result  in  an  error  value  of 

h?hr6^K,  where 


h “ time  step  sise. 


b ■ — for  Hamalng's  method. 


®^*i*  value,  h b 6 k,  was  small  %rtien  coapared  with  the  absolute  value  of 
the.  truncation  error,  l.e., 

< |t^|. 


Zn  all  the  computer  runs, 

h^b^a^K  - 0 , 


which  was  indeed  less  than  the  absolute  value  of  the  trwcatlon  error  in 
every  instance  of  the  calculations.  Also,  utilizing  this  procedure, 
only  two  iterations  of  the  corrector  were  required  for  convergence  in 

most  oases. 
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5.  BXPOHBMTXM.  DXmRBMCZMG  METHOD 

T)m  Moond  Mthod  ua«d  to  oolv*  difforontlal  c_A«tlons  (1)  and  (2) 
«aa  an  aaponantlal  dlffaranelng  tadmlqna  davalopad  by  Popa.^  Thla 
■athod  haa  baen  ahown  to  hava  superior  stability  propartlas  for  larga 
step  slsas  whan  daaling  with  a larga  class  of  dlffarantial  aquations. 
Hanca#  this  tachnlqua  asiy  ba  usad  with  a larqa  stap  sisa  to  decrease 
significantly  the  total  cmputing  tlsM  for  a aolutlon,  particularly  in 
those  anginaaring  problansr  liha  IMP  problMW»  whara  high  accuracy  is 
not  nacassarily  needed.  Bowavar,  in  this  work,  tlia  accuracy  is  a 
significant  part  of  tha  effort^  and»  as  it  turns  out*  the  axpcnantial 
differencing  aathod  does  prowida  precise  results,  whi^  are  shown  in 
section  7.  Tha  asponantial  differencing  technique  is  now  described. 

Oonsldar 

y'  - fUry)  O) 

with  tha  initial  condition  y(i^)  ■ y^.  The  ai^onantlal  diffaranoa 
equation  usad  to  solve  equation  (3)  is  given  by 

*2/  \ 

where 

h is  the  tiae  step  sise, 

y"  - f + f.f  , 
j X y 

f^,  fy  are  tha  partial  derivatives  of  f with  respect  to  x and  y.  Tha 
truncation  error  for  this  algorithm  is 

Vl”  ■ 

It  the  value  of  |hfy|  is  small,  at  least  if  |hf  | < 0.1,  than  in  place 
of  tha  exponential  formula,  the  serins  form  ^ 

- 

Vi  ■ if.  * “ * 'f”  ^ 


^Dsvid  A.  Pope,  An  Exponential  Method  of  Numerical  Integration  of 
Ordinary  Differential  Squatione,  Communications  of  the  ACm,  6,  No.  d' 
(Auguet  1963),  491-493. 
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is  ussd  to  avoid  loss  of  significance  due  to  cancellation  of  terms.  In 
this  case/  only  a few  terms  of  the  series  are  needed;  we  used  only 
three.  If  the  value  of  hf  is  fairly  large,  then  the  esqjonential 
subroutines  should  be  used.  ^ 

To  find  the  solutions  for  equations  (1)  and  (2)  using  the 
e^q>onential  differencing  method,  we  let 

^ - T ^ 

g(t,E)  ■ ^ ^ N^(t)E  , t > 0 . 

The  various  equations  needed  to  use  equations  (4)  and  (5)  are  now 
generated.  For  the  function  f(t,E),  its  derivatives  and  partial 
derivatives  are 

f _ - - 3^  H tOe"^  , 

E e e 

E"  - f^  ♦ f.fg  - f - quN^(0)E] 

+ N^C0)e2*^^N^(0)E  - J^]  , 

f^^  - e*^[j^  - <IMN,(0)e]  , 


f „ ■ 0 * 

EE 


f ^ M^(0)re**  . 


The  tnmeation  error  is 


‘n+l 


« i h>  fi  e*‘[j^  - q»«^(0)E]  +1^  «,(0)e2'^[qmi^(0)l  - jjj 
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For  th«  function  g(t,E) , its  dorivatluus  and  partial  darivatives  ara 
found  to  ba 

9 - - ft  »•  (t)  - - sat  s*E  , 

V C V C 

s--f  ".w  ' 


f 
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MOHERiaU.  VRLOBS 


Finding  thn  aolutlon  of  oqantionn  (1)  and  (2)  roqulrod  tho  on*  of 
tho  MM  hypothotieal  noMrieal  vnluoa  u i^liod  by  Wyntt^  for  tho 
various  oonstant  paroMtoro.  Amm  valuas  are  as  follows  i 

q - 1.6*10"*®  C, 

H^(0)  • 1.133*10“/oi3^ 

V m 2.0  m^/V/m  , 

J (t)  “ 1. 34*10^  for  all  t , 

c 

e • 8,85* 10“*  P/k  , 

X - 5.0* 10® /a  , 

B(0)  - 3. 697*  10**  V/m  . 

lha  fraa-alactxon  danslty«  (t) , is  givan  by 

Matt)  - 8^(0)  ♦ St  , (6) 

ahata  t ia  tha  tlM  and  8 is  tha  rata  of  production  of  ion  pairs  par 
unit  voluM  and  ia  glvan  tha  valua 

8 - 5.664*102®/k®/s  . 


Sinoa  aqoation  (6)  is  valid  for  aMthana  only  to  9 ns*  any  solutions  for 
tiMS  bayond  this  wara  not  posslbla  to  oaleulata.  Bosavar,  it  was 
poosibla  to  obtain  sKoallant  raaults  and  Mka  sqm  intarasting 
ooayarisons  for  this  xastrietiva  tlaa  frMa*  whi<di  ara  shown  ia 
saetion  7. 


r.  Atearaai  MF  atrmtrtk  ami  rim  nspsnisans  for  as 

fisaliaai  mtlrn,  J^poire  19M«  V.«.  Mhliitp  ffaifMae  UssssfcIi 

mi  asvalof—e  Osatar#  fbrt  Balmir,  m ^airaarp  1971), 
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RESULTS  AND  CX)MPARZS(»IS 


Tifo  conputer  codes  were  written  to  solve  differential  equations  (1) 
and  (2) . Coaiputer  code  HAMMING  was  written  to  solve  the  equations  by 
Hassling 's  method,  and  the  code  EXFDIFF  was  used  to  solve  the  equations 
by  the  eiqionential  differencing  Mthods.  The  listings  of  HAMMING  and 
EXPDXFF  are  in  appendices  A and  B,  respectively. 

Although  the  solutions  from  HAMMING  ooeqiare  favorably  «d.th  the 
results  of  Nyatt,^  the  solutions  obtained  are  considered  more  accurate, 
since  truncation  errors  are  calculated  and  examined.  in  table  I, 
specific  values  from  Nyatt^  and  HAMMING  are  listed  to  show  the  agreement 
of  the  results.  Regarding  these  results,  nuswrous  computer  runs  were 


TABLE  I.  RESULTS  CT  ORIGINAL  AMD  HAMIING 


Tiae 

(*) 

Original 
Electric  field 
strength 
(V/n) 

HAMMING 
Electric  field 
strength 

(At  - 1.0*10'”)  (V/b) 

-4.0.10"* 

3.80*10* 

3.6970*10* 

-J.0*10"’ 

3.80*10* 

3.6965*10* 

-2.0-10"* 

3.80*10* 

3.6961*10* 

-1.0* 10"’ 

3.80*10* 

3.6960*10* 

0 

3.80*10* 

3.6959*10* 

5.0* 10"'* 

3.00*10* 

3.2129*10* 

1.0* 10"’ 

2.48*10* 

2.6259*10* 

1.5*10"’ 

2.10*10* 

2.2092*10* 

2.0*10"’ 

1.85*10* 

1.9107*10* 

3.0*10"’ 

1.47*10* 

1.5093*10* 

4.0*10"’ 

1.20*10* 

1.2496*10* 

5.0* 10"’ 

1.00*10* 

1.0670*10* 

6.0* 10"* 

8.80*10* 

9.3132*10’ 

7.0*10"’ 

8.00*10’ 

8.2648*10’ 

•.0*10*’ 

7.40*10’ 

7.4296*10’ 

*.0*10"’ 

7.00*10’ 

6.7483*10’ 

r.  M^te#  Xmeamal  EXP  gtrmpgth  and  Tim  Dependenoe  for  m 
Xdaaiiaod  Trohlm,  Mmfort  IfM,  U,M,  Army  MoMlity  Egaijaant  Usasaiwh 
ami  Jmwaiofane  Caatar,  flart  Baltoir,  n (Fmbrmrr  1971), 


wmAm  with  dlffwrwnt  tia*  step  sisws,  which  ravwaled  sou  interesting 
facts.  When  the  tisw  step  size  At  - was  used,  stability  of 
the  solution  was  obtained  fron  0 to  approximately  6.5  ns,  tdiile  fr<»  6 
to  9 ns,  the  solution  became  dominated  by  very  large  truncation  errors 
and,  thus,  resulted  in  an  unstable  solution  over  this  time  frame.  Using 
the  tisM  step  sises  At  “ 1.0*10~^^  and  At  ■ 1.0*10~^^  did  not  produce 
any  stability  problesis  during  the  time  fram  0 to  9 ns,  and  ^ 
truncation  errors  were  considerably  smaller  than  for  At  - 1.0*10~^. 
The  truncation  errors  did  mt  decrease,  however,  tdien  the  time  step  size 
was  lowered  to  At  « 1.0*10  ^ , as  they  remained  approximately  on  the 
same  order  of  magnitude.  An  exas^le  of  these  results  appears  in 
table  II. 


TABLE  II.  HAMtlNG'S  TRUIKATION  ERRORS  FOR  DIFFERENT  TIME  STEP  SIZES 


Timm 

(•) 

Truncation  errora  for  tiae  step  sizes 

(At) 

t - i.e-io*'" 

Al  - 1.0*10*" 

t • 1.0*10*'' 

At  - 1.0*10*' ’ 

0 

0 

0 

0 

0 

1.0*  10'"’ 

0 

3.4909*10*’ 

-8.6590*10*" 

-1.2123*10"'* 

2.n*io''*' 

0 

4.1106*10*’ 

-1.7318*10*" 

-6.9272*10*' ' 

3.0-10'’' 

0 

4.0913*10*’ 

-5.1954*10"' ' 

-1.2123*10*'* 

4.0*10'"’ 

1.9735*10*' 

3.6824*10*’ 

-6.9272*10*" 

-1.3854*10*'* 

1.0- 10"’ 

8.5145*10** 

5.5282*10*’ 

-6.9272*10*' ' 

-6.0613*10*" 

2.0* 10'* 

-1.1649*10*’ 

-1.5846*10*’ 

-5.1954*10*" 

-4.3295  10*" 

3.0*10"* 

-1.7558* 10*' 

-3.2038*10*’ 

-6.4942*10*" 

-5.1954*10"" 

4.0-10*’ 

-1.0871*10*' 

-8.9188*10*'* 

-3.8965*10*' ' 

-2.5977*10*" 

5.0*10*’ 

-8.4018*10*’ 

-3.4203*10*'* 

-2.1647*10*" 

-3.0306*10*" 

6.0*10*’ 

-1.2692*10' 

-1.6019*10"'* 

-3.0306*10*" 

-3.0306*10"" 

7.0*10"’ 

-2.3989*10’ 

-8.6590*10*" 

-3.4636*10*" 

-3.0306*10"" 

8.0*10*’ 

-5.0072*10'* 

-6.4942*10*" 

-1.9483*10*" 

-2.1647*10"" 

9.0*10*’ 

-9.8303*10'’ 

-3.2471*10*" 

-1.2988*10*" 

-1.0824*10"" 

interesting  occurrences  were  noted  when  analysing  the  results 
of  nfOZrr.  First,  BXPDIFF  was  very  simple  to  code.  In  fact,  because 
of  the  relative  ease  and  simplicity  in  programing  BCPDIFP,  valid 
rmmmlts  ware  obtained  on  the  first  oompater  run.  Second,  the  results  of 
BFDXFF  were  quite  favorable  whan  oonpared  to  the  results  obtained  by 
iQfmtt*  and  again  are  oonsidered  mors  aooorate  because  of  the  calculation 


^If.  r.  Ifgaet,  XneemaJ  EKP  Strength  and  Tim  Dependence  tor  an 
Idealised  Problem,  Deport  1994,  0.8.  Army  Mobility  Equipment  Reeeareh 
and  Development  Center,  Port  Eelvoir,  VA  (Pebruarg  1971). 


TABLE  XZX.  RBSOLTS  OP  ORXGXNAL,  HAIMXIIG,  AND  EXPDXFP 


TIm 

(a) 

0 r 1 g i na 1 
Electric  field 
strength 
(V/m) 

HAHMING 
Electric  field 
strength 

At  - 1.0  • I0*‘2 
(V/n) 

EXPDIFF 
Electric  field 
strength 

At  - 1 .0  10*12 

(V/m) 

-4.0  • 10*’ 

3.80 

• 10‘ 

3.6970 

. 10* 

3.6970 

. 10* 

-3.0  • 10*’ 

3.80 

• 10‘ 

3.6965 

. 10* 

3.6965 

. 10* 

-2.0  • 10** 

3.80 

• 10‘ 

3.6961 

. 10* 

3.6961 

. 10* 

-1.0  . 10*’ 

3.80 

. 10’ 

3.6960 

. 10* 

3.6960 

. 10* 

0 

3.80 

. in'* 

3.6959 

. 10* 

3.6959 

. 10* 

5.0  • U'*” 

3.00 

• 10" 

3.2129 

• 10* 

3.3135 

. 10* 

1.0  • 10*’ 

3.48 

• 10" 

2.6259 

. 10* 

2.6764 

. 10* 

1.5  • 10*’ 

2.10 

. 10" 

2.209? 

. 10* 

2.2095 

. 10* 

2.0  • 10*’ 

1.85 

• lO’ 

1.9107 

. 10* 

1.9110 

. 10* 

3.0  • 10*’ 

1.47 

• 10- 

1.5093 

• 10" 

1.5094 

. 10* 

4.0  • 1C*’ 

1.20 

• lO*- 

1.2496 

. 10* 

1.2497 

. 10* 

5.0  • 10*’ 

1.00 

• 10* 

1.0670 

• 10* 

1.0670 

. 10* 

6.0  • 10*’ 

8.80 

• 10’ 

9.3112 

. 10’ 

9.3138 

. 10’ 

1 

O 

fH 

o 

8.00 

. 10’ 

8.2648 

. 10’ 

8.2652 

. 10’ 

8.0  • 10*’ 

7.40 

. 10’ 

7.4296 

. 10’ 

7.4299 

. 10* 

9.0  • 10*’ 

7.00 

• 10’ 

6.7483 

. in’ 

f.7486 

. 10’ 

•nd  analysis  of  tha  truncation  arrors  (tabla  XXX) . Probably  the  most 
signlfioant  aspect  of  BXPDXFF  ms  the  results  of  conpeurlsons  with 

lANMINa. 

Xa  oeaparlng  the  oodes,  several  significant  points  of  Interest  mre 
rsvaaled.  There  ms  a noticeable  ease  in  the  prograaalng  of  EXPDXFP 
ooNparsd  with  tha  aore  difficult  effort  required  for  HAMtlNG.  in  fact, 
a very  generous  analysis  of  this  coding  effort  ms  that  BANNING  took  at 
least  twice  as  long  to  code  as  BXPDIPP.  Althou^  HAIMXNG  is  generally 
oonsidared  to  give  acre  accurate  calculations,  EXPDXFP  yielded  results 
****  eNtreaely  close  (soawtiass  exact)  to  the  solutions  obtained 
froa  HANNING  (table  ZZZ).  However,  the  aost  significant  result  ms  the 
stabiUty  of  EXPDOT  when  ocaparod  with  HANNING.  For  tha  tiae  st^ 
sise  At  ■ 1.0*  10***,  tha  truncation  arrors  were  approalaately  the  saas 
else,  eaoapt  in  several  eases  where  truneatlmi  errors  froa  banning  were 
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•n  ordar  of  wignltudo  anllor.  Mian  tha  tlaw  step  alca  waa  inoraasad 
to  ^ HMMZMG  haraaa  unstabla  at  about  6.5  ns,  idiila 
BXPDIFF  yialdad  truncation  artors  on  tha  order  of  10~^  for  0 to  9 ns. 
This  tiadAg  is  lllustsatad  In  table  XV.  An  even  more  meaningful  result 
was  noticed  whan  oonsldarlng  tha  large  step  si^a  of  At  - 1.0*10~^.  For 
this  step  sisa,  BXPDIFF  was  stable  through  9 ns,  since  truncation  errors 
on  the  oirdar  of  10 wore  tolerated.  Whereas  BXPDIFF  exhibited 
reasonable  truncation  errors  for  this  step  slsa,  HAMMING,  as  expected, 
became  dominated  by  the  buildup  of  truncation  errors  and,  thus,  was 
unstable.  The  truncation  errors  of  BXPDIFF  for  At  ■ 1.0* 10~^  can  be 
seen  in  table  V. 


TABLE  XV.  TRDNCATXCM  BSBORS  kOR  DIFFERENT  TIME  STEP  SIZES 
FOR  HAMMING  AND  BXPDIFF 


Tloie 

(8) 

Truncation  errors  at  time  step 
size  1.0  • 10*‘“ 

Truncation  errors  at  time  step 
size  1.0  • 10*“ 

HAMMING 

EXPDIFF 

HAMMING  EXPDIFF 

0 

1.0  • 10'“ 

2.0  • lO"'* 

3.0  • 10'“ 

4.0  • lO"" 

0 

c 

0 

0 

1.9735  • 10*' 

7.4679  • 10*’ 
5.9718  * 10“" 
5.1673  • 10"* 
8.4422  • 10*‘ 
1.0322  • 10"* 

0 5.0168  • 10** 
3.4909  • 10**  3.9421  * 10"* 
4.1106  • 10**  6.740  • 10** 
4.0913  • 10**  8*3.3  * 10‘* 
3.6824  • 10**  9.0334  • 10"* 

1.0  • 10‘* 

8.5146  • 10*' 

9.7018  • 10** 

5.5282  *10*’  7.0506.10** 

2.0  • 10“* 

'1.1649  • 10*’ 

5.3903  • 10** 

-1.5846  • 10**  3.6655  * 10** 

3.0  • 10** 

•1.7558  • 10*‘ 

3.5640  • 10** 

-3.2038  • 10**  2.2718  • 10** 

4.0  • 10** 

-1.0871  • 10** 

2.6160  • 10** 

-8.9188  . 10*'*  1.5602  • 10** 

5.0  • 10** 

'8.4018  • 10*’ 

2.0477  • 10** 

'3.4203  r 10*"  1.1436  • 10** 

*.0  • 10** 

'1.2692  • 10 ' 

1.6749  • 10** 

'1.6019  • 10*"  8.7743  • 10*'* 

7,0  • 10** 

'2.3*89  • 10* 

1.4142  -10** 

'8.6590  • 10*"  6.9648  • 10*" 

8.0  • 10** 

'5.0072  • lo" 

1.2230  • 10** 

-6.4942  • 10*"  5.6760  • 10*" 

*.0  • 10** 

'9.8303  • lO'* 

1.0775  • 10** 

-3.2471  • 10*"  4.7242  • 10*'* 

TABLE  V. 


EXPDIFF  TRUNCATION  ERRORS  AT  TIME  STEP  SIZE  1.0  ’ 10' 


-9 


Time  (s) 


Truncation  errors 


0 

8.7131 

10"' 

1.0 

10"’ 

2.1944 

10"’ 

2.0 

10”’ 

5.1673 

10"’ 

3.0 

10"’ 

3.4598 

10"* 

4.0 

10"’ 

2.5851 

in'* 

5.0 

10*’ 

2.0672 

10'* 

6.0 

10”’ 

1.7226 

10"’ 

7.0 

10'* 

1.4766 

10"* 

8.0 

10"’ 

1.2920 

10"* 

9.0 

10"’ 

1.1484 

10‘* 

8.  CONCLUSIONS* 

With  regard  to  the  caaparlsons  made  between  the  predictor-corrector 
routine  HAMMING  and  the  ejgponential  differencing  method  EXPDIFF,  there 
sees  to  be  two  cooqpelling  factors  that  m2Jce  the  exponential  differencing 
technique  8iq>erior.  First,  the  ease  and  simplicity  of  the  prograasdng 
effort  required  for  EXPDIFF  far  outweigh  the  more  coiqplicated  coding 

work  needed  for  HAMMING.  Second  and  probably  more  isqportant,  the 
stability  properties  of  EXPDIFF  are  excellent,  tdiereas  the  solutions 
calculated  by  HAMMING  became  dominated  by  truncation  errors  during  the 
fraste  examined  for  particular  time  step  sizes.  This  stability 
property  is  best  exeaqplified  by  the  large  step  size  (1.0*10~^)  that  can 
be  used  with  relative  assuredness  of  accurate  results.  Thus,  this 
stability  factor  r^resented  by  the  reasonable  triincation  errors  of 
EXPDIFF  far  outweigh  the  someidtat  smaller  truncation  errors  of  HAMMING. 
Therefore,  considering  the  ease  of  programming  and  the  stability 
properties,  the  exponential  differencing  method  provides  an  efficient 
and  reliable  solution  to  the  differential  equations  involved  in  this  EMP 
problem  and  is  strongly  recommended  for  solving  other  differential 
equations  of  this  type. 
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